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Abstract. We propose damage spreading (DS) as a tool to investigate the 
topological features related to the ruggedness of the free energy landscape. 
We argue that DS measures the positiveness of the largest Lyapunov exponent 
associated to the basins of attraction visited by the system during its dynamical 
evolution. We discuss recent results obtained in the framework of mode-coupling 
theory and comment how to extend them to the study of realistic glasses. 
Preliminary results are presented for purely repulsive soft-sphere glasses. 
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Introduction 

The theoretical understanding of a first principles theory for the glass transition is 
still missing. Despite of great advances in the understanding of some generic features 
associated to the glass transition (such as those predicted by the mode-coupling 
theory) still some questions remain largely unknown. Going beyond the schematic 
mode-coupling theory seems to be an enormous task so an alternative way of looking 
at the glass transition may be useful. In this direction, the study of the topological 
properties of both the potential or free energy landscape may yield further information 
on the mechanisms responsible for the anomalous viscosity of the glassy phase. 

The idea that topological aspects of the potential or free energy landscape are the 
ultimate reason for the glass transition goes back to Goldstein Q and more recently 
Stillinger and Weber ^ ||. This approach has been recently applied to the study of 
hard-spheres Q] , monoatomic as well as binary Lennard- Jones glasses [|| or mean- field 
models of glasses || . 

Here we propose an alternative dynamical approach to study the topological 
properties of the potential energy landscape. We will concentrate on the study of 
the stability local properties of the configurations visited by the system during its 
dynamical evolution. This is directly achieved through the study of how dynamical 
trajectories, which evolve following the same stochastic noise, depart from each 
other in the presence of a potential energy saddle point or a maximum which may 
induce a negative Lyapunov exponent. The simplest way to study this problem is 
through damage spreading (DS) techniques to be describe later on in some detail. 
Although DS was introduced almost two decades ago as an alternative way to consider 
thermodynamic phase transitions, the initial enthusiasm on this problem strongly 
decayed when it was realized that DS transitions are not universal and not necessarily 
related to thermodynamic singularities. 

Despite of this result here we will show that these transitions have an added 
interest in that they may be used as a direct way to investigate the local free energy 
landscape properties by measuring the largest Lyapunov exponent associated with the 
Hamming distance (to be defined later). In what follows I will explain in more detail 
why DS is a good way of looking at the rugged properties of the potential energy 
landscape. Later on I will discuss the analytical results obtained for the schematic 
mode-coupling theory and finally discuss how to extend these ideas to the study of 
real glasses. Some preliminary results are shown for the case of binary soft-sphere 
purely repulsive glasses. 

1. Why damage spreading? 

Consider two systems evolving under a Langevin dynamics each one described by a 
set of N variables Xi,yf,l < i < N evolving in a potential energy landscape V under 
the same stochastic noise rj t with (r]i(t)r]j(s)) — 2T5ijS(t — s). Although the present 
discussion can be generalized for different stochastic noises here we will concentrate 
on the simplest case (for a more detailed discussion see Q). The equations of motion 
read, 



x t (t) = Fi({x}) + rn(t) 
yi(t)=F i ({y}) + r H (t) 



(1) 
(2) 
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where Fi({x}) = Note that both trajectories described by the systems x and y 

never cross in phase space so two identical configurations such that x t (t) — (t) remain 
identical forever (and where identical in the past). The equation for the difference 
variables z, = Xi — yi reads, 

Zi(t) = Fi({ Xi }) - Fittyi}) . (3) 

If the Zi are small we can expand (j^) around Zi = obtaining, 

which may be written in a simplified form, 

Zi = Hij({y})zj (5) 

where Hij is the Hessian matrix evaluated at the configuration y. Always within the 
linear approximation the dynamical evolution of the distance between configurations 
Zi will increase or decrease according whether the spectrum of eigenvalues of the 
Hessian matrix contains positive eigenvalues. In this sense, DS probes the spectrum 
of eigenvalues of the matrix and shows instabilities whenever the matrix develops 
positive eigenvalues. A more precise condition is given by the maximum Lyapunov 
exponent defined through, 

X max = lim (6) 

t — >oo t 

where D{t) = jt zf which should be positive whenever Zi — is dynamically 
unstable. Note that the Hessian depends on time through the time evolution of 
the generic configuration y. This may be an equilibrium or an off-equilibrium 
configuration. So in principle the maximum Lyapunov exponent depends on time 
through the time evolution of the systems x and y. We will see later that, in general, 
the type of initial condition (as well as the initial distance) are not relevant parameters 
for the DS transition. In this sense DS probes the temperature at which the lowest 
accessible configurations in the potential energy landscape develop unstable modes 
being a direct check of the corrugated properties of the free energy landscape. Again, 
we must insist on the non-universal properties of the DS dynamics. The present 
discussion on the stability properties of the Hessian matrix and its connection with 
the DS transition is valid in the framework of Langevin dynamics. For other type of 
dynamics (such as Monte Carlo dynamics or Glauber) the situation may be different 
and the physical meaning of DS phenomena more difficult. In some sense, Langevin 
dynamics is an appropriate tool to explore the topological properties of the potential 
energy landscape. 



2. DS in mode-coupling theory 

Insight on the previous problem can be obtained through a careful study of the DS 
equations in the case of ideal mode-coupling theory. It is known since the seminal work 
by Kirkpatrick, Thirumalai and Wolyness || that mode coupling equations can be 
obtained in the framework of exactly solvable p-spin glass models. Due to their mean- 
field character, in this class of models it is possible to unambiguously define concepts 
such as the configurational entropy or complexity and the mode-coupling transition 
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temperature T c . The description of this type of models is possible in the framework 
of the TAP analysis || where it is possible to show that they contain a large number 
of metastable states (exponentially large with N) as well as a threshold energy where 
the system gets trapped in an asymptotic aging state and the fluctuation-dissipation 
theorem is violated in a peculiar way jig ]. 

Spherical p-spin models (contrarily to Ising spins) have the clear advantage of 
being exactly solvable so it is convenient to do analytical computations in that case. 
The potential energy in this model is defined by, 

V = - X] Ji 1 ,i2,i 3 ,..,i p < J ti°'i2< J i3--< J ip ■ (7) 

(il<i2<---<i P ) 

where the spins &i are real valued spins which satisfy the spherical constraint 
2i=i °f = N . The Ji 1 ,i 2 ,i3,--,ip are quenched random variables with zero mean and 
variance p!/(2A rp_1 ). The Langevin dynamics of the model is given by, 

^l=F i ({a})- m + r ]i (8) 

where [i is a Lagrange multiplier which ensures that the spherical constraint is 
satisfied at all times and the noise r\ satisfies the fluctuation-dissipation relation 
(Vi(t)Vj( s )) — 2T5(t — s)5ij where (...) denotes the noise average. Fi is the force 
acting on the spin Oi due to the interaction with the rest of the spins, 

Fl = = („ _ l)\ X] Jii,i2,--,i P (T i2 Cr i3-- a i P ■ (9) 



We define the overlap between two configurations of the spins a, r by the relation 

N 



Q = jf Ylj—i °~iTi so a distance between these two configurations is 



D=^ (10) 

in such a way that identical configurations have zero distance and opposite 
configurations have maximal distance D = 1. Then we consider two copies of the 
system {o~i,Ti} which evolve under equation (^|) with the same statistical noise and 
start from random initial configurations. 
The final equations are , 

' K "" s ' 1 fi(t)C(t, s) = | / duR(s,u)C p - 1 (t,u) + 



dt r-v v v w 2 j o 



p(p - 1) 

2 



duR{t,u)C(s,u)C p - 2 (t,u) (11) 







^ifl + fi(t)R(t, s) = S(t - s) + ^P—^L f duR(t,u)R(u,s)C p - 2 (t,u) (12) 

Ot Z J s 

+ ^(t)Q(t, s) = \f Q duR(s, u)Q p -\t, u) 

+ p(p_pl) f duR (t )U )Q(u lS )C p - 2 {t,u) (13) 
2 Jo 

The dynamical equations involve the two times correlation, response and overlap 
function C{t, s),R(t, s), Q(t, s) defined by (in what follows we take t > s), 
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N N 

C(t, s ) = (i/tf) $><(«)*•<(*)) = (i/N)J2(n(t)Ti(s)) (u) 
i=i »=i 

R(t, S )=(i/N)f:^=(i/N)f:^ as) 

z—l 1 i—1 % 

N 

Q(t,s) = (l/N)J2(^(t)n(s)) (16) 

i=l 

where < .. > denotes the average over dynamical histories and hf,hj are 
fields coupled to the spins <7j,T, respectively. These equations are complemented 
with the appropriate boundary conditions C(t,t) = l,Q^(t) = Q(t,t), R(s,t) = 
0, linit— »( s )+ R(t,s) — 1 and the relations, 



i i(t) = T + ^- [ duR(t,u)C p - x {t,u) (17) 



ldQd ^ . tl (t)Q d {t)=T+l [ duR{t,u)Q p -\t,u) 



2 at • 2 

ft 



(I 



P{j> ^ f duR(t,u)Q(t,u)C p - 2 (t,u) . (18) 



These equations can be analyzed in detail using different methods. Here we 
summarize the main results obtained Q], 



Existence of a dynamical transition To 

There is a temperature To such that D{t) = (or Qd(t) = 1, see eq.(|l0|)) is a stable 
fixed point for T > Tq becoming unstable below To. Because of the non-monotonic 
character of D(t) it is very difficult to derive analytically Tq. Nevertheless, it is 
possible to obtain an upper and a lower bound. One gets, 



^ < To < (19) 

Direct numerical integration of the equations of motion yields T (p = 3) = 
1.04 ± 0.02 with and T (p = 4) = 1.13 ± 0.02. The value of T is well above 
the mode-coupling temperature T c and the TAP temperature Ttap below which 
there is an exponentially large (with the system size) number of metastable states. 

Independence of initial conditions 

The asymptotic damage D(po) — lim^oo limjv^oo D(t) is independent on the 
value of the initial damage D(0) or the class of initial conditions (for instance, 
random or thermalized). This independence stresses the fact that DS is a true 
dynamical transition and the asymptotic damage D(oo) is a dynamical order 
parameter. 

Tq is the lowest DS temperature 

The DS problem can be suitably generalized for the case of correlated noises such 
that (r]i(t)£,j( s )) — 2TJC(Q(t, s))6(t — s)Sij where r\ and £ are the noises acting 
on the systems a and r respectively. The function /C satisfies JC(1) = 1 so both 
noises are identical if the two configurations coincide. This implies that Qd(t) = 1 



Dynamical phase transitions in glasses induced by the ruggedness of the free energy landscaped) 




Figure 1. Asymptotic distance Doc for p = 3 (a = 1, K = 1) obtained from 
the Pade analysis of the series expansions for different initial conditions Do = 1 
(circles), Do = 0.5 (triangles), Do = 0.25 (stars). Typical error bars are shown 
for the last case. 




Figure 2. Qd(t) for p = 3 (a = 1,IC = 1) at temperatures T = 0.1,0.5 (from 
bottom to top at large times) for three different values of the initial overlap 
Qd(0) = —1,0,0.5 as a function of time. The continuous lines are the numerical 
integrations with time step At = 0.01. 



is a fixed point of the dynamics. It can be shown that for any possible function 
K < 1 (with JC(1) = 1) there is a finite temperature damage spreading transition 
To only if /C'(l) < 1. The case discussed previously K, — 1 (identical noises at all 
times) yields the lowest damage spreading transition temperature. 

• T is the endpoint of a dynamical critical line 
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The DS problem can be also generalized to the case fC(Q) < A with A < 1 and 
/C(l) = A. Obviously for A = 1 dynamical trajectories of both systems may 
cross. In this case it is possible to show that the function JC(Q) = A yields 
the lowest DS transition temperature To (A) among the set of possible functions K, 
{K.(Q) < A, /C(l) = A). T (A) is monotonic increasing function function of A which 
for A = coincides with the mode-coupling transition temperature T c and finishes 
in a critical endpoint To (A = 1) = Tq. So there exists a line of dynamic critical 
points which connect the mode-coupling temperature T c with the DS temperature 
T . 

• T) is not universal. 

The temperature To is not universal. As it depends on the set of correlations of the 
noises it also depends on the type of dynamics (molecular dynamics, Monte Carlo 
with Metropolis, heat-bath or Glauber). This is a well known result which finds 
its natural explanation on the physical origin of the DS transition. For a general 
dynamics it is not possible to map the DS transition with the local properties 
of the potential energy landscape. Only for the case of Langevin dynamics 
or molecular dynamics this is possible. Other dynamics (such as Monte Carlo 
with heat-bath dynamics) use random numbers in the dynamics which introduce 
complex correlations between the noises. This yields a DS transition (related to 
the To (A) discussed in the previous paragraph for the Langevin case) which is 
probably related with the mode-coupling transition temperature but this issue 
still needs to be further investigated. 



3. Application to binary soft-sphere glasses 

In this section we apply the previous ideas derived in the framework of mode-coupling 
theory to the case of structural glasses. We consider the binary soft-spheres model 
introduced in and recently studied in [|l2| . For sake of simplicity we consider a 
gas of N particles such that half of them have diameter o\ and the other half 02. The 
particles interact through a two particle purely repulsive potential, the energy of the 
system being defined by 




The choice ct^ = a ' ^ supposes that diameters are additive during the collision 
process. The advantage of this potential is that the thermodynamic properties 
depend on the density p = N/V and the temperature T only through the constant 
T = p/Ti. For the particular case ^ = 1.2 crystallization is strongly inhibited and 
the glass transition (where dynamics is strongly slowed down) appears in the vicinity 
of r = 1.45. Larger values of V correspond to the glass phase while lower values 
correspond to the liquid phase. The Langevin dynamics for the soft-sphere model is 
defined by, 

N 

= -^2^ l V lJ (r ll )+rf l (21) 

with {rii{t)ri l j{t')) = 2TSijSkiS(t — t') where the superindex in the noise indicate the 
different Cartesian components of the vector noise ff(t). The pairwise potential is 
given by V %3 (r) = (^) 12 . 
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We now consider two systems described by the variables fl , s\ governed by ( pl[ ) 
and evolving under the same realization of the noise. We define the Euclidean distance 



which vanishes if the two configurations coincide. If we want to extend the 
previous ideas for the spherical p-spin model to this system now we must take into 
account the fact that at very high temperatures a gas diffuses so D = may not be a 
fixed point of the dynamics. There are two strategies to deal with this problem which 
are discussed below. 

• Particles contained in a box 

This is the most natural choice. To simulate a purely repulsive system one must 
confine the particles in a cubic box of side L such that p = pp In this case 
one may numerically solve ( |2(i| ) with two different class of boundary conditions. 
With periodic boundary conditions particles leave one side of the box and enter 
the opposite side. This resets completely the coordinates of the particle so the 
distance (p2| ) is discontinuous if particles cross the boundaries. Concerning one 
system quantities (such as the energy or the pair correlation function) this is not 
a problem because the relevant quantity is the distance between the particles 
which may be taken as the minimum value between ry and L — rij. A similar 
procedure can be used to define the distance between the two copies. Everything 
can be easily solved considering free boundary conditions so particles are not 
allowed to cross the boundaries. In this case, it is possible to show that D = is 
asymptotically stable for the purely diffusive case (T = 0). 

Preliminary results show that the DS transition temperature Tq — oo so two 
configurations never coincide at finite temperature. Still both configurations 
retain some correlation (so (fi(t)sl(t)) > 0) and the asymptotic damage is a 
non trivial function of the temperature. 

• Introducing an spherical constraint 

For the purpose of studying the local properties of the potential energy landscape 
we may impose the global constraint J2i ^ — ^(ir) s 011 the particles in such 
a way that the average distance between the particles is finite when N goes 
to infinity. Because the spherical constraint shifts the Hessian matrix (|^) by a 
constant (a Lagrange multiplier) the transition with the spherical constraint may 
give information on the transition for the unconstrained case. That Lagrange 
multiplier can be simply obtained from the potential energy < V > and the 
temperature. The advantage of such a constraint is that now there is no box and 
D = is a fixed point of the dynamics for T = 0. The inconvenience is that the 
simplicity of the original model is lost and the thermodynamics of the new model 
depends on both density and temperature instead of a unique parameter T. 
Again, preliminary results show that To = oo in this case so D = is 
asymptotically stable strictly only for T = 0. Although this approach is more 
involved it is probably the best way to relate the DS transition to the ruggedness 
of the free energy landscape. 




(22) 



Dynamical phase transitions in glasses induced by the ruggedness of the free energy landscaped 




Figure 3. Damage D(t) as a function of time for N = 32 starting from two 
different initial conditions and three different temperatures (from top to bottom) 
T = 0.8,0.4,0.2 



4. Conclusions 

The study of the free energy landscape may yield valuable information on the glass 
transition phenomena. A promising description of the glass transition is through 
the Stillinger and Weber projection of the partition function in terms on inherent 
structures. That method directly looks at the potential energy landscape described in 
terms of basins of attraction explored by the system during its dynamical evolution 
||. An alternative approach studies the dynamical properties of the free energy 
landscape directly looking at the largest Lyapunov exponent of the Hessian matrix 
of the potential energy landscape weighted by the size of the basins of attraction 
visited by the system during its dynamical evolution. 

Exact results for the mode coupling theory reveal that there is a transition To 
which separates two well defined regime depending on the value of the asymptotic 
distance. Below To the asymptotic damage is non zero and independent of the initial 
distance as well as the class of initial conditions. Above To the damage vanishes. We 
argue that the precise value of T is related to the vanishing of the largest Lyapunov 
exponent defined in (||). Although such an explicit connection needs still to be done 
it is quite probable that DS is a precise tool to investigate the chaotic properties 
of the free energy landscape. A result in this direction has been recently obtained 
by Biroli through the study of the instantaneous normal modes spectra of the p- 
spin model Whether this transition has experimental relevance in the study of 
real glasses is still an open question. Our preliminary studies of soft-sphere binary 
mixtures show that Tq is extremely large. Because liquids are always diffusive at large 
temperatures (a feature which is directly encoded in the wave-vector dependence of 
correlation functions, a general feature of liquids) one must be careful when extending 
the results obtained for the spherical p-spin model to real structural models of glasses. 
Although a better understanding of the extension of DS to diffusive systems is needed 
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we can point out other interesting open problems. One the one hand it could be 
very interesting to analyze the DS transition for molecular dynamics. In that case, 
there is no stochasticity in the dynamical equations so the effective source of noise 
comes out directly from the mixing property of the dynamics. The analog of equation 
should be very similar except for the presence of oscillations. Still the general 
argument would be the same and To expected to be identical. Such an analysis would 
be welcome. Finally it would be very interesting to look at the other endpoint of the 
dynamic critical line. Our present discussion was centered on the case of identical 
noises. For completely uncorrelatcd noises the dynamical transition temperature is 
expected to coincide with the mode-coupling transition temperature. This is true in 
the framework of the aforementioned exact calculations in the spherical model and 
could be also analyzed for real glasses. 
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